PROCEEDINGS 

OF SCIENCE 



Phase Transition Signature Results from PHENIX 



J.T. Mitchell^ 

Brookhaven National Laboratory 
E-mail: }nitchell@bnl . gov 



The PHENIX experiment has conducted searches for the QCD critical point with measurements of 
multiplicity fluctuations, transverse momentum fluctuations, event-by-event kaon-to-pion ratios, 
elliptic flow, and correlations. Measurements have been made in several collision systems as a 
function of centrality and transverse momentum. The results do not show significant evidence 
of critical behavior in the collision systems and energies studied, although several interesting 
features are discussed. 
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1. Introduction 

Recent work with lattice gauge theory simulations indicate that the phase diagram of Quantum 
Chromodynamics (QCD) may contain a first-order transition line between the hadron gas phase 
and the strongly-coupled Quark-Gluon Plasma (sQGP) phase that terminates at a critical point [jl]]. 
This property is observed in many common liquids, including water. Near the QCD critical point, 
several thermodynamic properties of the system will diverge with a power law behavior in the 
variable e = (T — Tq) /Tc, where Tc is the critical temperature. Here, several measurements made 
by the PHENIX experiment at Brookhaven National Laboratory's Relativistic Heavy Ion Collider 
that may be sensitive to this critical behavior are discussed. 

2. Multiplicity Fluctuations 

In the Grand Canonical Ensemble, the variance and the mean of the particle number, N, can 
be directly related to the compressibility, kj: (On = var ^' = k B Tykj, where kg is Boltzmann's 
constant, T is the temperature, and V is the volume [§]. Near the critical point, the compressibility 
diverges with a power law behavior with exponent y: kj <=c s~ r . The measurement of event-by- 
event fluctuations in the multiplicity of charged hadrons may be sensitive to critical behavior in the 
system. PHENIX has surveyed the behavior of inclusive charged particle multiplicity fluctuations 
as a function of centrality and transverse momentum in y / ^v=62.4 GeV and 200 GeV Au+Au 
collisions, and in ^A AW =22.5, 62.4, and 200 GeV Cu+Cu collisions. 

Since multiplicity fluctuations are well described by Negative Binomial Distributions (NBD) 
in both elementary [^] and heavy ion collisions [Q], the data for a given centrality and pj bin are 
fit to an NBD from which the mean and variance are determined. Due to the finite width of each 
centrality bin, there is a non-dynamic component of the observed fluctuations that is present due 
to fluctuations in the impact parameter within a centrality bin. The magnitude of this component 
is estimated using the HIJING event generator which well reproduces the mean multiplicity of 
RHIC collisions The estimate is performed by comparing fluctuations from simulated events 
with a fixed impact parameter to events with a range of impact parameters covering the width 
of each centrality bin, as determined from Glauber model simulations. The data are corrected to 
remove the impact parameter fluctuation component. 

Baseline comparisons are made to the participant superposition model, in which the total mul- 
tiplicity fluctuations can be expressed in terms of the scaled variance [01, (On = (O v + Hwn <%»a»> 
where (O v are the fluctuations from each individual source, (% „ are the fluctuations of the number 
of sources, and jU^jv i s tne mean multiplicity per wounded nucleon. The second term includes non- 
dynamic contributions from impact parameter fluctuations along with additional fluctuations in the 
number of participants for a fixed impact parameter. Ideally, the second term is nearly nullified after 
applying the previously described corrections, so the resulting fluctuations are independent of cen- 
trality as well as collision species. Baseline comparisons at 200 GeV are facilitated by PHENIX 
measurements of charged particle multiplicity fluctuations in minimum bias 200 GeV p+p colli- 
sions with mean pL = 0.32 ± 0.003, scaled variance (O = 1.17 ± 0.01, and NBD fit parameter k^BD 
= 1.88 ±0.01. 
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Figure 1: Multiplicity fluctuations as a function of N par t for Au+Au collisions for 0.2 < pj < 2.0 GeV/c. 
Contributions from impact parameter fluctuations have been removed. Shaded regions represent a 1 a range 
of the superposition model prediction derived from p+p data. 



The scaled variance as a function of the number of participating nucleons, N part , over the pj 
range 0.2 < pj < 2.0 GeV/c is shown in Figure p] for Au+Au collisions. For all centralities, the 
scaled variance values consistently lie above the Poisson distribution value of 1.0. In all collision 
systems, the minimum scaled variance occurs in the most central collisions and then begins to 
increase as the centrality decreases. A similar centrality-dependent trend of the scaled variance has 
also been observed at the SPS in low energy Pb+Pb collisions at ^/smm=11.3 GeV, measured by 
experiment NA49 [JSj] , where the hard scattering contribution is expected to be small. All of the 
data points are consistent with or below the participant superposition model estimate. This suggests 
that the data do not show any indications of the presence of a critical point, where the fluctuations 
are expected to be much larger than the participant superposition model expectation. 

The clan model []|] has been developed to interpret the fact that Negative Binomial Distribu- 
tions describe charged hadron multiplicity distributions in elementary and heavy ion collisions. In 
this model, hadron production is modeled as independent emission of a number of hadronic clus- 
ters, N c , each with a mean number of hadrons, n c . The independent emission is described by a 
Poisson distribution with an average cluster, or clan, multiplicity of N c . After the clusters are emit- 
ted, they fragment into the final state hadrons. The measured value of the mean multiplicity, jj^, is 
related to the cluster multiplicities by /i c h = N c n c . In this model, the cluster multiplicity parameters 
can be simply related to the NBD parameters of the measured multiplicity distribution as follows: 



N c = kmo log(l + Hch /kmr>) (2.1) 

and 

n c = (jUchANBD)/fog(l +aW&nbd)- (2-2) 
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Figure 2: The correlation of the clan model parameters n c and N c for all of the collision species measured 
as a function of centrality. Also shown are results from pseudorapidity-dependent studies from elementary 
collisions (UA5 [§], EMC PL and NA22 [O]) and heavy ion collisions (E802 [1 and NA35 pi]). 



The results from the NBD fits to the data are plotted in Fig. |2] for all collision species. Also 
shown are data from elementary and heavy ion collisions at various collision energies. The indi- 
vidual data points from all but the PHENIX data are taken from multiplicity distributions measured 
over varying ranges of pseudorapidity, while the PHENIX data are taken as a function of centrality. 
The characteristics of all of the heavy ion data sets are the same. The value of n c varies little within 
the range 1.0-1.1. The heavy ion data universally exhibit only weak clustering characteristics as 
interpreted by the clan model. There is also no significant variation seen with collision energy. 
However, n c is consistently significantly higher in elementary collisions. In elementary collisions, 
it is less probable to produce events with a high multiplicity, which can reveal rare sources of 
clusters such as jet production or multiple parton interactions. 



3. (pr) Fluctuations 

PHENIX has also completed a survey that expands upon previous measurements of event- 



by-event transverse momentum fluctuations Q13Q. Here, the magnitude of the pj fluctuations will 



be quoted using the variable L PT , as described in []15|]. H PT is the mean of the covariance of all 
particle pairs in an event, normalized by the inclusive mean pj. Z pr is related to the inverse of the 
heat capacity of the system [jl6|], which diverges with a power law behavior near the critical point: 

Cy «£-«. 

Figure |5| shows L PT as a function of N part for all 5 collision systems measured over the pj range 
0.2 < pj < 2.0 GeV/c. The data is shown within the effective PHENIX azimuthal acceptance of 
4.24 radians. The magnitude of E PT exhibits little variation for the different collision energies and 
does not scale with the jet cross section at different energies, hence hard processes are not the 
primary contributor to the observed fluctuations. Simulations show that elliptic flow contributes 



little [13]. With the exception of the most peripheral collisions, all systems exhibit a universal 
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power law scaling as a function of N part . The data points for all systems are best described by 
the curve: Z P7 oc N pc } r ^ 2±0A0 . The observed scaling is independent of the pj range over which the 
measurement is made. 
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Figure 3: Event-by-event pj fluctuations for inclusive charged hadrons within the PHENIX acceptance in 
the transverse momentum range 0.2 < pr < 2.0 GeV/c in terms of T. PT as a function of N par t- 



4. Kin Fluctuations 

PHENIX has studied identified particle fluctuations by measuring the event-by-event fluctua- 
tions of kaons to pions and protons to pions. One advantage of particle ratio measurements is that 
contributions from volume fluctuations cancel. Measurements are quoted in the variable Vd yn '■ 



(n (K-l)) (K(K-\)) (Kit) 
<7T>2 (Kf (K)(n) 



v dyn (K,n) = + (4.1) 



If only random fluctuations are present, Vd yn is zero. Also, Vd yn is independent of acceptance. 

The measurements for Vd yn (K,7l) for 0.34 < pj < 1.05 GeV/c are shown in Figure |j. The 
measurements for Vd yn (K,p) are shown in Figure |5| As with the pj fluctuations, the fluctuations 
in (K) / (n) demonstrate a VN part dependence. This is not seen in fluctuations of (p}/(n), which 
instead rise as centrality increases. 

5. Scaling of Elliptic Flow 

One of the most striking RHIC results has been the observation of scaling behavior in elliptic 
flow measurements below a transverse momentum of about 1 GeV that indicate that quark degrees 
of freedom are driving the dynamics of the collision Jl4|]. PHENIX measurements of the scaling 
behavior of elliptic flow are compiled for various particle species and various collision systems in 
Figure |6| Further measurements of this scaling behavior and the observation of its breaking as a 
function of collision energy will be an important ingredient in the search for a critical point. 
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Figure 4: Event-by-event fluctuations of the kaon-to-pion ratio for inclusive charged hadrons within the 
PHENIX acceptance in the transverse momentum range 0.35 < pr < 1 05 GeV/c. The dashed line is a fit to 
the function c+Np^, where c is a constant. 
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Figure 5: Event-by-event fluctuations of the kaon-to-proton ratio for inclusive charged hadrons within the 
PHENIX acceptance in the transverse momentum range 0.35 < pr < 105 GeV/c. 



6. Searching for a Critical Point with HBT Correlations 

Near the critical point, correlation functions are also expected to be described by a power law 
function with critical exponent tj. This exponent can be measured with Hanbury-Brown Twiss 
correlations in the Qi m , variable [|l8|]. Here, the Q\ m correlations are fit with a Levy function, 

C{Q im ) = Xexp{-\Rq/hc\- a ), (6.1) 

where R is the HBT radius and a is the Levy index of stability. The value of a is 1 for a Lorentzian 
source and 2 for a Gaussian source. Since a equates to the exponent tj, it is expected that its 
value will approach the value expected for the universality class of QCD. If QCD belongs to the 
3d Ising model class, the value of tj would approach 0.5. Figure [7] shows the results of the Levy 
function fit to PHENIX HBT correlations in 0-5% central 200 GeV Au+Au collisions as a function 
of transverse mass. The fit results are inconsistent with the expected value of 0.5 in the vicinity of 
a critical point. Analysis of the other PHENIX datasets is currently underway. 



7. Azimuthal Correlations at Low Transverse Momentum 

Critical behavior may also be apparent in the width and shape of correlation functions. PHENIX 
has measured azimuthal correlation functions of like-sign pairs at low pj for several collision sys- 
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Figure 6: PHENIX Preliminary elliptic flow V2 normalized by the number of quarks, collision eccentricity, 

1 /3 

and N part plotted as a function of the transverse kinetic energy normalized by the number of quarks. Shown 
are V2 measurements of pions, kaons, and protons over a centrality range of 0-50% for 200 GeV Au+Au, 
200 GeV Cu+Cu, and 62.4 GeV Au+Au. 
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Figure 7: The Levy index of stability a extracted from Levy function fits to g™> correlations in 0-5% central 
200 GeV Au+Au collisions as a function of transverse mass. 



terns. The correlations isolate the HBT peak in pseudorapidity by restricting \Ar] \ < 0.1 for each 
particle pair. Correlations are constructed for low pj pairs by correlating all particle pairs in an 
event where both particles lie within the pj range 0.2 < pj \ < 0.4 GeV/c and 0.2 < pj2 < 0.4 
GeV/c. Note that there is no trigger particle in this analysis. The correlation functions are con- 
structed using mixed events as follows: C(A<f>) = ™l*Z daa N '~^ . Confirmation of the HBT 
peak has been made by observing its disappearance in unlike-sign pair correlations and by observ- 
ing Qinvariant peaks when selecting this region. 

Azimuthal correlation functions can be described by a power law function with exponent tj : 
C(A0) oc A(j>~( d ~ 2+T1 \ where d is the dimensionality of the system For all collision systems, 
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including 200 GeV d+Au, the extracted value of the exponent 77 is shown in Fig. ||. The value 
of v\ lies between -0.6 and -0.7 with d=3, independent of centrality. Since v\ is constant in heavy 
ion collisions, does not differ from the d+Au system, and has a value that significantly differs from 
expectations from a QCD phase transition (e.g. tj=+0.5 for the 3-D Ising model universality class 
[17]), it is unlikely that critical behavior is being observed in the correlation functions measured 
thus far. 

Near the critical point, it is also expected that the correlation length will diverge with a power 
law behavior. The HBT peak of the correlation functions with the estimated contribution from 
elliptic flow subtracted have been fit to a Gaussian distribution. The standard deviation from the fit 
is shown in Figure |9] for several collision species. There is no significant change in the correlation 
widths between 200 GeV Au+Au and 62.4 GeV Au+Au collisions. 
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Figure 8: The exponent 77 with d=3 extracted from the like-sign correlation functions as a function of N part . 



8. Conclusions 

The fluctuation and correlation measures presented here do not provide a significant indication 
of the existence of a critical point or phase transition. This does not rule out the possibility that 
the critical point exists. Further searches will be facilitated by the upcoming RHIC low energy 
program. 
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